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The application of the exact renormalisation group to a many-fermion system with a short- 
range attractive force is studied. We assume a simple ansatz for the effective action with 
effective bosons, describing pairing effects and derive a set of approximate flow equations for 
the effective coupling including boson and fermionic fluctuations. The phase transition to a 
phase with broken symmetry is found at a critical value of the running scale. The mean-field 
results are recovered if boson-loop effects are omitted. The calculations with two different 
forms of the regulator was shown to lead to a similar results. 
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There is a growing interest in applying the Exact Renormalisation Group (ERG) formalism 


to few- and many-body systems 


L apDlymg_i 


when the underlying interaction is essentially non- 
perturbative. Regardless of the details all ERG-based approaches share the same distinctive fea¬ 
ture, a successive elimination/suppression of some modes, resulting in effective interaction between 


the remaining degrees of freedom j^. One specific way of implementing such a procedure is to 
eliminate modes by applying a momentum-space blocking transformation with some physically 
motivated cutoff. The effect of varying a cutoff is described by nonlinear ERG evolution equations, 
which include the effect of the elliminated modes. By solving the ERG equations one can find a 
scale dependence of the coupling constants and thus determine a path in the space of Lagrangean 
functionals. The ERG formalism is closely linked to an other approach which has become increas¬ 
ingly popular in both few- and many-body physics, Effective Field Theory (EFT). EFT is also 
based on a separation of scales, removal of some (mainly high-energy) degrees of freedom and use 
the effective degrees of freedom instead of the fundamental ones. 

In a sense, EFT and ERG compliment each other. EFT can provide the guidance for fixing 
the initial conditions and ERG can be used as an analytic method to study the evolution of the 
system as a function of some generic scale parameter. This is especially useful in the case of truly 
non-perturbative problems, where there are no small parameters one can expand in. This situation 
is quite common both in few- and many-body problem. One notes that, although we will focus on 
the systems consisting of nucleons, many aspects of the following discussion are relevant for the 
other types of fermionic systems, especially fermionic atoms in traps. 
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Probably the most important dynamical featnre of the nucleon-nucleon interaction is the un¬ 
naturally large scattering length that makes a perturbative expansion meaningless both for free 
nucleons and in nuclear matter. In addition, the presence of the Fermi momentum signals the 
appearance of another scale which further complicates the use of any perturbative technique. One 
notes that the use of the perturbation theory for the two-nucleon system in vacuum can be justi¬ 
fied only for the hypothetical case of weakly interacting nucleons with a small scattering length.In 
medium, however, even a weak attraction between nucleons may lead to the intrinsically nonper- 
turbative phenomenon, the superfluidity, characterised by rearrengement of the ground state and 
appearance of the gap in the spectrum. The fermions form correlared pairs which, depending on the 
strength of the interaction, may lead to different physical regimes. The weak coupling regime (BCS 
phase) corresponds to a pair with the spatial size much larger then the radius of the interaction 
so that no actual bound two-body subsystem is formed, while in the strong regime corresponding 
to the Bose-Einstein Condensation (BEC) the fermion pairs form compact deeply bound two-body 
states. 

It would be an appealing to describe all these regimes by starting at a large initial scale with 
some EET motivated effective action and then run the scale down to the physical point where the 
scale parameter vanishes. Thus the results will depend on one physical parameter, the scattering 
length of the nucleon-nucleon interaction in free space. 

The ideal tool to treat this problem is provided by a variant of the ERG approach, based on 
the average effective action (AEA) 

irreducible (IPI) correlation function in the presence of an infrared cutoff scale k. Only fluctuations 
with momenta larger then k are taken into account. Eor /c —> 0 all fluctuations are included and we 
arrive at standard effective action from which all physical correlation functions can be extracted. 
The evolution equation for the AEA has the following one-loop form 
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Tr 


{dkR) (r(2) - R) 



( 1 ) 


Here E^^) is the second functional derivative of the AEA taken with respect to all types of field 
included in the action, and i? is a regulator which should suppress the contributions of states with 
momenta less than or of the order of running scale k. To recover the full effective action we require 
R{k) to vanish as /c —> 0, in other respects its form is rather arbitrary. The concrete functional 
form of the regulator has no effect on physical results provided no approximations/truncations were 
made. In practice, however, approximations/truncations are always required to render the system 
of the evolution equations finite and solvable. Therefore, some dependence on the functional form 
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of regulator is inevitable. The simplest way to estimate this dependence is to solve the system of 
evolution equations using several choices of regulators. Of course, it does not garantee the fully 
quantitative estimate of the errors introduced but may at least give an idea about the size of the 
corresponding uncertainty. 


I. ANSATZ FOR T 


We demand that at high scale our theory be a purely fermionic theory with the contact inter^ 
action decsribed by the lagrangian 
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Since we are interested in the appearance of the correlated fermion pairs in a physical ground state, 
we need to parametrise our effective action in a way that can describe the qualitative change in 
the physics when this occurs. A natural way to do this is to introduce a boson field whose vacuum 
expectation value (VEV) describes this correlated pair |3|.At the start of the RG evolution, the 
boson field is not dynamical and is introduced through a Hubbard-Stratonovich transformation of 
the four-point interaction. As we integrate out more and more of the fermion degrees of freedom 
by running k to lower values, we generate dynamical terms in the bosonic effective action. In this 
paper we treat a single species of fermion. The corresponding ansatz for the boson-fermion effective 
action can be written as 




(fx (j)\x) ( {idt + 2n) + ^ (j){x) 
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( 3 ) 


Here M is the mass of the fermions in vacuum and the factor l/2m in the boson kinetic term 
is chosen simply to make dimensionless. The couplings, the chemical potential /r, the wave- 
function renormalisations and the kinetic-mass renormalisations Zm,M all run with on k, the 
scale of the regulator. 

The bosons are , in principle coupled to the chemical potential via a quadratic term in cp, but 
this can be absorbed into the potential by defining U = U — 2p.Z^(pP(p. 

We expand this potential about its minimum, cp'^cp = po, so that the coefficients Ui are dehned 
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at /) = /5o, 

11 

U{p) = uo + ui{p - po) + ■7:U2ip - pof + -^U‘i{p - pof H-, (4) 

1 o 

where we have introduced p = The phase of the system is determined by the coefficient ui. 
In the symmetric phase we have po = 0 so that the expansion takes the form 

U{p) =Uo + Uip+]^U2P^ ^ -. ( 5 ) 

The potential in the condensed phase can be simplified to 

U{p) = uo + ]^U2{p-pof+ ■■■ . ( 6 ) 

In our current work we shall truncate this potential at quartic order in the field (order p^). However, 
the fact that we define our coupling constants at the minimum of the potential does mean that we 
need to consider the next term in the expansion. This will allow us to treat the implicit dependence 
of the coefficients on pQ. We treat the wave function renormalisation factor for the bosons in the 
same way, expanding it about p = pQ as 

Z^{p) = z^Q +z^i{p-Po) ^ -. (7) 

The other couplings and renormalisation factors can be treated similarly. 

The fermions are not dressed at this point and the bosons are just auxiliary fields and so we 
can assume that Z^{K) = 1, Zm{K) = 1. 

We have a choice between following the evolution for fixed chemical potential, and allowing 

p to run so that the fermion density is kept fixed. The region corresponding to a Bose-Einstein 

condensate (BEC) of tightly bound pairs corresponds to negative values of p and is not accessible 
to evolution at fixed chemical potential. We therefore follow the evolution at fixed density, n. In 
this context it is convenient to define the Eermi momentum, pp, corresponding to this density and 
the chemical potential at the starting scale by 

PF = p{K) = (8) 

This connection between p and pp only holds in the symmetric phase, where p does not run. 

II. EVOLUTION EQUATIONS: GENERAL STRUCTURE 

In this section we consider the general structure of the evolution equations. We start by con¬ 
sidering the simpler case of evolution at constant chemical potential. The boson potential U is 
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obtained by evaluating the effective action for uniform boson fields. It is given by 


dkU = -^dkr, 

1/4 


(9) 


where V 4 is the volume of spacetime. Substituting our expansion of U, Eq. on the left-hand 
side leads to a set of ordinary differential equations for the Un- We choose to evolve following the 
minimum of the potential and use the expansion around the minimum p = po{k) to define the Un, 
as in Eq. ©• In the symmetric phase, p = 0, it gives the following set of equations 

dun 


dk (9p" 


dkU 


p=0 


In the condensed phase we need to define the total derivative 

— = f) + ^ — 
dk ^ dk dpo 


( 10 ) 


( 11 ) 


Acting on d'^U/dp^ with this (and taking the higher derivative term over to the LHS) gives a set 
of equations, each of which is coupled to the coefficient of the next term through the evolution of 
P 01 


dur 




dpQ 5 " 




( 12 ) 


P=P0 


dk dk dp'"' 

The coefficient ui is special since it vanishes for the expansion around the minimum. Imposing the 
condition ui{k) = 0 gives an equation for the evolution of po; 


«!&(/ 

dk op 


(13) 


P=P0 


The coefficient uq provides information on the energy density of matter. It satisfies the equation 

(14) 


^^0 _ o rn 
dk 


which does not couple back into the other equations. 

In the condensed phase, we could truncate our potential at quadratic order, as we do on the 
right-hand sides of these equations, and simply set U 3 = 0 on the left hand-side of the equation 
for U 2 - However, a better approximation can be obtained by substituting the exact form for u^{k) 
taken from the evolution with fermion loops only, as described in below. This is the approach we 
adopt here. It has the benefit of providing an approximation to U 2 {k) and po{k) that becomes 
exact in situations where boson loops can be neglected. 

The evolution of the boson wave-function renormalisation factor can be obtained from 

1 9 / 52 \ 

V4 dpo 


( 15 ) 
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If we substitute our expansion, Eq. 0 , we get another set of coupled equations in the condensed 
phase. Only the first of these is of interest within our current truncation 


dz^o _ ^ 
dk dk 


J_ A 

V4 dpo 




po=0, p=po 


(16) 


Again, corresponds to a term beyond our current level of truncation and so we will take the 
result from fermion loops only. The evolution of can also be deduced from 


dkZtj, 


1 

2 dpdp 



(17) 


which gives 


dz^o _ ^ 
dk dk 
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2 dpdp 



P=P0 


(18) 


The evolution equations for the other couplings {Zm, Zm, Zp, Zg) can be derived in a similar 
manner. However, in this paper we allow to run only Z^, parameters in the potential and chemical 
potential since this is the minimal set needed to include the effective boson dynamics and study 
the BCS-BEC crossover. 

The fermion number density is given by 


n = — 


dU 


dp 


P=P0 


The evolution equation for n can be written as 


(19) 


The equations constructed so far describe the evolution at constant p. If we want to follow 
the evolution at constant density we must allow p to run with k. In this case we define the total 
derivative 


( 20 ) 


A - ^ _L d 

dk ^ dk dpo dk dp' 

Applying this to dtJ/dp at p = po gives the evolution equation for n 

dn „ dpn dp d 




P=P0 


where we have introduced the fermion-number susceptibility 


d‘^U 

dp^ 


( 21 ) 


( 22 ) 


P=P0 


(23) 
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If n is kept constant {dn/dk = 0) this becomes 


dpo dn _ d / - 

^ dk^^dk~ dp 


(24) 


P=P0 


This equation describes evolution of chemical potential in the broken phase 

The remaining set of evolution equations to be solved (in the broken phase) is 


duo , dp 

—rr+n— = OkUl _ , 

dk dk ' P-Po 


dpo 

dk 


-U2 — + 2z<^o 


dU2 

~dk 


dpo 

dk 


U3 -JJ- + 2Z(j,i 


dp 

Ik 

dp 

lik 


where we have dehned 

, dW 

X = 


dp^dp 


A 

dp 

dp'^ 




dkU 


P=P0 


P=P0 


dz^o dpo I , dp 
~dk 2^ Im 


dzmO dpo ^ dp 



? 7=0 


p2=o, p=pQ 


d^ ( 52 


OLrm. — 


dpT 


P=P0 


V4 dpd{p^) \dpdp^ ^ 


p=0 


p2=0, p=po 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


The set of evolution equations in symmetric phase can easily be recoved using the fact that chemical 
potential does not run in symmetric phase and that po = 0. The left-hand sides of these equations 
contain a number of coefficients that lie beyond our current level of truncation, such as x, U 3 and 
We propose to replace these by their exact expressions obtained from evolution with fermion 
loops only. The formal derivation of the corresponding expressions will be considered in more 
details below. 


III. CHOICE OF CUT-OFF 


In the bosonic sector, we take the regulator to be an additional quadratic term, proportional 
to (j)^ {x)(j){x'). In the representation used to write down the second derivatives above, it has the 
matrix structure 

/ RB(q,k) 0 \ A:2 

^B{q, k) = , Rp,{q,k) = -^ f{q/k), (31) 

V 0 RB{q,k) ) 2m 

where f{x) ^ 1 as x ^ 0 . 
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In the fermion case, our regulator should be positive for particle states (Zm9^/2M > Z^fi) 
and negative for hole states {ZmQ^/‘^M < Z^p). It should suppress the contributions of states 
with energies near One easy way to ensure this would be to use the off-diagonal regulator and 
so generates an artificial gap in the fermion spectrum around fi. However such a regulator could 
not be used without a Fermi sea and so would not allow us to connect our results in matter to 
the interaction between the fermions in vacuum. We therefore choose our regulator to have the 
structure 


RF{q,PF,k) = 


where we have introduced 


/ 


sgn{q-p^)RFiq,PF,k) 0 

0 -sgn{q - Pf,) RF{q,PF,k) 


= 


I Z^2Mp 
Zm ’ 


(32) 


(33) 


the Fermi momentum corresponding to the (running) value of p. 

The function RF{q,PF,k) should suppress the contributions of states with momenta near the 
Fermi surface, \q — pf\ ~ k. Once a a large gap has appeared in the fermion spectrum, there are 
no low-energy fermion excitations and so the fermionic regulator plays little further role. However, 
while the gap is zero or small, it is crucial that the sign of the regulator match that of the energy, 
ZMq‘^/‘^M — Z-^p, and hence it is p which appears in the sign functions. 

In order for to be sure that we are matching onto the same bare NN interaction at the starting 
scale K, the fermionic regulator should satisfy RFiq,PF, K) — RF{q,0, K)- 


IV. DRIVING TERMS: POTENTIAL 


In this section we derive the evolution equations for the parameters of the potential. Calculating 
the fermion propagator, multiplying by B^Rf and taking the matrix trace gives 


where 



{dkRF) (rg]. 



2EFR{q,PF,k) sgn{q - p^)dkRF{q,PF,k) 
^p(lo - EFR{q,PF, kf - A‘^ + ie 


(34) 


^M 

EFR{q,PF, k) = — q^ - Z^p -h RF{q,PF, k) sgn(g - p^), 


The poles in this propagator occur at 


qo = E—y^EFR{q,PF, ky + A2. 


(35) 


( 36 ) 
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At A: = 0 (iip = 0) in the condensed phase, these become 

and so the gap in the fermion spectrum at q = Pf is 2A/Z^. 

The corresponding boson matrix trace can be worked out in a similar way. Putting everything 
together gives 


dkU = -^dkT 

1/4 


+ 


1 


(fq 

(27r)3 

f 

(2f)^ 


Efr 


^FR + 

Err 


sgn{q - p^) dkRF 


^BR 


dkRB- 


(38) 


where 


EBR{q,k) 



2m 


q'^ + ui+ tt2(2(/)V - po) + RB{q, k), 


Vb = n2</>V- 


(39) 


The driving terms in the evolution equations for the coefficients in our expansion, Eq. 0 , are 
obtained from the derivatives of dkU with respect to p = (p^cj) so that we get in the condensed 
phase 


A 

dp 



P=P0 


2Z. 


(fq 

(27r)3 


Efr 


EIr + A(^)2 


sgn(g - Pf,) dkRF 


+ 


U2V^^ 

2Z^ 


d^q 

(2f)^ 


Ei,i-2VP 

i 7'('^)2 T/(c)2^ 

Err ^B 


dkRB, 


(40) 


where 


E^^BRi^l) = ^<1^+ '^'iPQ +RB{q,k)]Vlf’=U2Po, Ai'=i=p^, (41) 

and similarly for the second derivative of dktJ with respect to p = (/)^(/). The equations for the 
couplings tti( 2 ) in the both phases can be obtained in a similar way from the driving terms remem¬ 
bering that in the symmetric phase (p = po,rii 7 ^ 0). The driving term for the evolution of the 
fermion number density is given by the derivative of U with respect to p. The sign functions in the 
fermion part of Eq. depend on p and so, in principle, differentiating with respect to p could 
generate surface terms. However the sign change occurs at precisely the point where Err vanishes 
and hence such terms do not arise. The resulting driving term vanishes in the symmetric phase, 
and so evolution at constant n and at constant p are the same there. 
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For the evolution of fj, in the condensed phase we get 




sgn(g - Pp) dkRp 




dkRs- 


V. DRIVING TERMS: Zj, 


For the wave function renormalisation factor, Z^, we need to consider a time-dependent 
ground field. The evolution of is then driven by 

1 d / I 


Dehning 


V4 dpo ydpdrj'^ 


■p(3) _ ^ ■p(2) _ 

J- DTDJ. - , -L DD - 


lpo=0 

—2u2(j)^ —2u2(t> 
0 —2u2(I)^ 


r(3) _ t.(2) _ I 0 ^ 9(^2 


we can write the relevant part of the evolution equation in the form 


dpdrj^^ 


dkV = +i Tr [(4Rf) (rgj, - Rf)-^ , (rgg - Rf)-^ rgg. (rgg - Rf)- 


^ Tr [{dkKB) (rgg - Rf)-' rggg (rgg - Rf)"' rg% (rgg - Rf)- 


After lengthly algebra one can obtain in broken phase 

1 d f f (fq 2El^ - A(^)2 


f Sq 24t-(54°>,'4°> + FW" 


77 '(c )2 -rr(c)2 

^BR 


dkRs- 


The corresponding expression in symmetric phase follows rather trvially. 



11 


As a check on this result, we note that ui contains a piece —2^Z^. Hence we can also obtain 
the evolution of from 

1 92 


2 d^dp 


U 


(48) 


P=PO 


Taking the partial derivative with respect to p, as discussed at the end of the previous section, the 
result we obtain agrees with Eq. ( 1171 ) . 

VI. FERMIONS ONLY 


In order to estimate the effects of the higher order coefficients (us etc) we use the ERG equations 
when boson loops are neglected. In this case the expressions simplify considerably and the total 
effective potential can be calculated analytically. All the needed coefficients can then be expracted 
by simple differentiation. The RG equation for the effective potential becomes 

d^q Efr 


dkU = - 


(27r) 


^FR + 


sgn{q-pf,) dkRp- 


(49) 


(In the absence of boson fluctations, we have Zp = 1.) This can be rewritten in the form 

d?q 


dkU = -di 


/ 


(27r) 


^FR. + > 


and immediately integrated to give 

d^q 


Uip,p,k) = Uip,p,K)- 


(27r)^ 


y^EFR{q,PF, kY + - y/EFR{q,PF, 


(50) 


(51) 


We have made explicit the dependence of the potential on the chemical potential p. 
At our starting scale K, we take the potential to have the form 

U{p,p,K)=uo{K) + ui{K)p. 


(52) 


The renormalised value of tti(i^)can be deduced from the scattering length. To determine uq{K), 
we use the fact that at p = 0 the physical potential is just that of a free Fermi gas, measured 
relative to the chemical potential: 

d^q 


and hence 


U{0,p,0) = 2 


f Q f dj^ Q 

uo{K) = J EFR{q,PF,0) - J j^;^EFRiq,PF,K)sgn{q-pF) 


(2vr) 


■ EFR{q,PF,0)e{pF - q), 


d^q 


(53) 


( 54 ) 
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The physical potential (at fc = 0) is then given by 

cfq 


Uip,n,0) = 


(27r)3 


EFRiq,PF, 0) — \/ EFRiq,PF, 0)^ + 


MA2 A2 f d^q 1 

47ra 2 


(55) 


(27r)3 £;^^(gf,0,0) ' 

Differentiating with respect to p and setting the derivative equal to zero, we find that A^ at the 
minimum satisfies 


+ 



1 

I 

1 (2vr)3 

EFR{q,o,o) 

\/EFR{q-,PF,^y + 


= 0 . 


(56) 


This is exactly the gap equation used in Ref.j 

To get the number density of fermions, we can differentiate U{p,p,0) with respect to p. This 
gives 


n = 


[ d^q 

^ EFR{q,PF,0) 

I (27r)3 

y/EFR{q,PF,0)^ + 


(57) 


in agreement with Ref. 


VII. INITIAL CONDITIONS 

We start our evolution at some large, fixed value for the cut-off scale, K, and evolve down to 
k = 0. This clearly requires some initial conditions on the parameters in our effective action. Some 
constraints on these have already been mentioned in Sec.Q but we now turn to fixing specific values 
for them. These should be derived from the known interactions between the particles in vacuum. 
In the vacuum case, it is relatively straightforward to fix the initial values for the parameters at 
the starting scale K. 

In matter, the detemination of our renormalised parameters is complicated by the fact that our 
fermionic cut-off depends on the Fermi momentum. The values of parameters such as ui{pf,K) 
must thus depend on pp as well as K. Note that we assume that our starting scale K is sufficiently 
large that all physical effects of the Fermi sea have been completely suppressed by our cut-off. Hence 
the dependence of the initial parameters on pF is merely to compensate for the pF dependence 
of our regularisation and renormalisation procedures. For example, in the case of a simple sharp 
cut-off, the maximum momentum included is K-\-pF in matter instead of K. It is thus particularly 
important to correctly renormalise ui{pF, K), since this quantity cancels a linear divergence and 
so can be shifted by a finite amount even for K —> oo. 
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One could introduce a cut-off function that tends to a p^-independent form ioi K ^ pp. How¬ 
ever in practice a modification of the renormalisation procedure is more convenient. In the region 
K ^ pf, we can ignore boson loops. The evolution of quantities such as ui{pf,K), U 2 {pf,K), 
Zff,{pF, K) and Zm{pF, K) is thus similar to the vacuum case, except for the different cut-off. This 
allow us to define ui{pF-, K) to be 

ui{pF,K) _ M l r d^q 

g‘^ Apa 2 J {2p)^ 

This expression can be thought of as being generated by the vacuum evolution using a modified 
cut-off that interpolates smoothly between Rpiq^PF-, k) for k ^ pp and Rf{q, 0, k) for k < pp. It 
ensures that our renormalised parameter ui{pf,K), defined using Rp^q^pp^k) for large k, corre¬ 
sponds to the physical scattering length in vacuum. 

The initial values for U 2 {pf,K), Zfj,{pF,K) and Zm{pF,K) can be determined using similar 
procedures, although this is not so crucial since these quantities do not contain linearly divergent 
pieces and so all their piT’-dependence is suppressed by powers oi pf/K. One convenient choice is 
to take their starting values to be zero at some large but finite scale K. An alternative is to require 
that they tend to zero as K ^ oo. 

The initial condition for the energy density is most conveniently expressed in terms of ho which, 
in the symmetric phase, is simply given by the energy of a free Fermi gas, measured relative to the 
chemical potential, and so its initial value is just 

uo{K)=2 J -^^EFR{q,PF,Q)0{pF - q)- (59) 

VIII. RESULTS 


I _ sgn(g - pp) 

I .Ei7’_R(g, 0,0) EFR{q,PF.,K)\ 


(58) 


We solve the evolution equations numerically with two types of cutoffs. First, we use the 

smoothed step-function type of regulator (called hereafter as Ri) 

^2 ^2 

Rif = ^0i{q - pF,k,a), Rib = ^0i{q,k,a) (60) 

where 

1 


ei{q,k,a) = 


2Erf(l/cj) 


Erf(i±f)+Erf(?-''' 


(61) 


ka ka 

with a being a parameter determining the sharpness of the step. 

Second, we use a sharp cutoff function denoted R 2 chosen to make the calculations as simple 
as possible 

R 2 F = ^ [((A: + p^,f - q‘^)9{Pi, + k-q) + {{k+ p^f + q^ - 2p^)6l(g -p^, + k)] , 


(62) 
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FIG. 1: Numerical solutions to the evolution equations for infinite oq and pf = 1-37 fm, starting from 
K = 16 fm“^. We show the evolution of all relevant parameters for the cases of fermion loops only (thin 
lines), and of bosonic loops with a running (thick lines). All quantities are expressed in appropriate 
powers of fm“^. 


R 2 B = - q). 


(63) 


Similar boson regulator was used in Ref. Q] (see also Ref. ^). 

As we can see the fermion sharp cutoff consists of two terms which result in modification of 
the particle and hole propagators respectively. The hole term is further modified to suppress the 
contribution from the surface terms, which may bring in the dangerous dependence of the regulator 
on the cutoff scale even at the vanishingly small k. As an example, we focus on the parameters 
relevant to neutron matter: M = 4:.76fm~^ ,pF = 1.37/m“^. 

We first discuss the results obtained with the smooth cutoff Ri. We find that the value of 
the physical gap is practically independent of either the values of the width parameter a (varied 
within some range) or the starting scale K provided K > The results of the calculations 

are shown on Fig. 1. At starting scale the system is in the symmetric phase and remains in this 
phase until ui hits zero at kcrit — 1.2fm~^ where the artificial second order phase transition to 
a broken phase occurs and the energy gap is formed. Already at A; ~ 0.5 the running scale has 
essentually no effect on the gap. We found very small (on the level of 1%) contribution to the gap 
from the boson loops, due to cancellations between the direct contributions to the running of the 
gap and indirect ones via U 2 - The boson loops play much more important role in the evolution 
of U 2 and Z^. In fact, they drive both couplings to zero at A; —> 0 although at rather slow pace. 
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FIG. 2: Evolution of the chemical potential 

We note however, that the effect of the boson loops for the gap may still be more visible if the 
evolution of the other couplings is included. The results obtained for the gap correspond to the 
case of the infinite negative scattering length. To study the BCS-BEC crossover we have to solve 
the evolution equation for a wider range of the scattering lengths, including the positive values. 
The corresponding results for the evolution of chemical potential as a function of the parameter 
Pfcl are shown in Fig.2. While the vacuum scattering length is large and negative, the system is 
in the BCS-like phase with positive chemical potential, whereas if the scattering length is chosen 
to be large and positive reflecting the existence of a bound state near threshold the system ends 
up being the collection of weakly overlaping tightly bound pairs with negative chemical potential. 

Now we turn to the results obtained with sharp cutoff (Fig.3). One immediate observation 
is that that the results become starting scale independent as long as K > 5fm~^ similar to the 
smooth regulator case. However, the phase transition occurs at lower values of the running scale 
k ~ 0.7/m“^. At approximately k ~ 0.2/m“^ the value of the gap becomes scale independent. 
The gap evolutions obtained with the smooth and sharp regulators, being rather different at inter¬ 
mediate scales, approach each other with decreasing scale resulting in similar values for the physical 
gap. This is an encouraging results taking into account that, although the exact results must be 
independent of the choice of the regulator, in practice it is not garanteed. The same conclusion also 
holds for other quantities. The couplings and U 2 first grow with scale and then start decreasing 
eventually coming to zero. Chemical potential begins to decrease at the point of phase transition 
and becomes scale independent at A: ~ 0.2fm~^. However, in this case the numerical values of 
the chemical potentials obtained with different regulators differ by approximately 20% so that 
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FIG. 3: Evolution of the parameters when the sharp cut-off is used 


this quantity is more sensitive to the details of effective action and to the trancations made. We 
note that the sharp regulator can also describe the BCS-BEC crossover, although gives somewhat 
smaller (negative) values of chemical potential at ppa > 1. Applying our approach to neutron 
matter we hnd a gap com par able to ep, of the order of 30 MeV. There is a simple explanation 
for the smaller values (see |l2l| , Q] or ) 
coupling, where the gap satisfies 


The argument can be given most succinctly for weak 


A = (8/e^)eFexp (- (7r/2) cot(5(pF))) • 


(64) 


For nucleon-nucleon scattering, cot 6 increases relatively quickly with momentum and the resulting 
reduction in the gap is substantial. We therefore expect that an extension of our approach to 
include the effective range should capture this physics. Indeed, if the “in-medium” scattering length 
is identified with the Bethe-Goldstone G matrix calculated at zero energy but finite momenta [1^ 
then we obtain the gap ~ 8 MeV which is already compatible with the commonly accepted value. 
Of course, this is only a crude estimate and the proper calculations should be done using the 
momentum dependent part of the four-fermion interaction. 

There are several points where our approach could be improved. We have already mentioned 
above the effective range effects. We should also include running of all the couplings and treat 
explicitly the particle-hole channels (RPA phonons) since these contain important physics. They 
will allow us to remove the “Fierz ambiguity” associated with our bosonisation of the underlying 
contact interaction We would like to include the three-body force effects, which are required 
to satisfy the reparametrisation invariance theorem and possibly the long-range forces. As 







17 


to further applications we plan to explore the superfluidity of the cold fermionic atoms in traps, 
temperature dependence of the BEC-BSC transitions and the formal relations of the ERG with 
the other many-body approaches. 

It is a pleasure to thank Mike Birse, Niels Walet and Judith McGovern for very useful discus¬ 
sions. 
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